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ABSTRACT 

This  report  has  been  produced  i  n  order  to  address  aerodynamic  characteristics  and  stabi  lity  and 
control  issues  relating  to  hypersonic  vehicles  that  are  deemed  to  be  significantly  different  from 
those  of  the  conventional  (subsonic  or  supersonic)  airvehides.  In  particular  we  have  addressed 
issues  that  are  relevant  to  stability  and  control  of  hypersonic  vehicles.  This  report  should  add  to 
the  existing  knowledge  of  missile  guidance  and  control  engineers  and  make,  other  researchers 
and  engineers  involved  in  hypersonic  experimentation,  aware  that  these  vehicles  may  not  be 
dynamically  stable  and  require  active  control  augmentation  in  order  to  achieve  and  maintain 
d  esi  rabl  e  f  I  i  ght  character  i  sti  cs. 
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On  Stability  and  Control  of  Hypersonic  Vehicles 


Executi  ve  Su  mmary 

Hypersonic  flight  presents  major  challenges  to  airframe  and  control  system  designers. 
High  velocity  can  cause  a  hypersonic  vehicle  to  be  highly  sensitive  to  changes  in  flight 
conditions  that  can  result  in  instability  or  weakly  damped  transient  oscillations  of  the 
airframe.  The  design  problem  is  further  compounded  by  the  fact  that  hypersonic 
aerodynamic  parameters,  as  predicted  from  ground  tests  or  theoretical  computational 
meth  od  s,  d  o  n  ot  ref  I  ect  th  e  actu  al  f  I  i  ght  p  arameters;  th  ere  are  si  gn  i  f  i  cant  u  n  certai  nti  es  i  n 
the  parameter  values  required  for  airframe  and  control  system  design.  Conventional 
techniques  do  not  always  lead  to  a  design  that  is  stable  and  at  the  same  time  robust  to 
parameter  uncertainties.  From  the  reported  work  in  hypersonic  vehicle  stability  and 
control,  we  are  able  to  highlight  aerodynamic  characteristics  that  may  have  a  strong 
bearing  on  our  approach  to  the  analysis,  synthesis  and  performance  evaluation  of 
hypersonic  weapons  systems.  The  key  findings  and  recommendations  are  as  follows: 

•  Active  control  will  be  required  to  maintain  stability  in  case  of  changes  to  the  CG 
position  and  maintain  adequate  damping  during  flight. 

•  Uncertainties  in  aerodynamic  parameters  require  a  control  system  design  based  on 
robust  techniques.  Adaptive  control  techniques  may  be  required  to  maintain  desirable 
vehicle  flight  performance. 

•  At  hy  personi  c  speeds  the  d  rag  and  I  ift  forces  become  non-l  i  near  fu  ncti  ons  of  the  angl  e 
of  attack.  Compared  to  subsonic  and  supersonic  speeds,  the  maximum  value  of  the 
I ifty  drag  ratio  for  hypersonic  vehicles  is  significantly  lower. 

•  H  eating  has  a  negative  affect  on  vehicle  structural  integrity  and,  depending  on  the 
airframe,  may  cause  structural  vibration  which  needs  to  be  catered  for  by  an 
appropriate  control  system  design. 
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Nomenclature 


M#  ; 

M  ach  number 

a  : 

Angle  of  attack 

y: 

Flight  path  angle 

9  -a  +  y : 

Total  body  angle 

CP  : 

Centre  of  pressure 

CG: 

Centre  of  gravity 

L: 

Aerodynamic  lift  force 

D: 

Aerodynamic  drag  force 

T : 

Thrust  force 

M : 

Aerodynamic  (pitching)  moment 

V : 

Vehicle  velocity  along  the  flight  path 

W : 

Vehideweight 

Cm  • 

Moment  Coefficient 

* 

q  : 

Dynamic  pressure 

S  : 

Characteristic  aerodynamic  area 

c : 

Characteristic  aerodynamic  chord  length 

p: 

Air  density 

CL: 

Li  ft  coefficient 

XCp  ■ 

D  i  stance  of  Cp  from  the  nose 

X(  (,  • 

D  i  stance  of  CG  from  the  nose 

XSM  ‘ 

Static  margin 

PC. 

Derivative  of  the  lift  coefficient 

c  =  L 

^ La  - 

da 

C  _  dCM 

D eri  vati  ve  of  the  1  i ft  coeff i  ci  ent 

M  a  ^ 

da 

T  : 

Xw 

xw  -  component  of  the  thrust  vector 

T  : 

zw 

zw  -  component  of  the  thrust  vector 

m  : 

Vehicle  mass 

8 : 

Earth's  gravity 

?».  •• 

<L  • 
q: 

I  : 

)7 

(Ow,xw,zK): 

(Ov,xv,zv): 

(variable)  : 

3 (variable) : 

1: 

//: 


Flight  path  rotation  rate  about  the  ytv  -axes 

Earth's  rotation  projected  along  the  yw  -axes 

Vehicle  pitch  rate  about  the  body  y  -axes 

Vehicle  body  moment  of  inertia  about  its  j-axes 

x,  z  coord i  nates  of  the  wi  nd-axes  system 

x,  z  coordinates  vehicle-axes  system 

Valueof  a  variable  at  a  given  operating  condition 

Deviation  of  a  variable  about  a  given  operating  condition 

Earth  latitude 
Earth  longitude 
Vehicle  control  deflection 
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1.  Introduction 

Hypersonic  flight  presents  major  challenges  to  airframe  and  control  system  designers.  High 
velocity  can  cause  a  hypersonic  vehicle  to  be  highly  sensitive  to  changes  in  flight  conditions 
(Mach  Number  M#,  and  angle  of  attack  a)  that  can  result  in  instability  or  weakly  damped 
transi  ent  osci  1 1  ati  ons  of  the  ai  rf  rame.  The  desi  gn  probl  em  i  s  fu  rther  compou  nded  by  the  fact 
that  hypersonic  aerodynamic  parameters,  as  predicted  from  ground  tests  or  theoretical 
computational  methods,  do  not  reflect  the  actual  flight  parameters;  there  are  significant 
uncertainties  in  the  parameter  values  required  for  airframe  and  control  system  design. 
Consequently,  conventional  techniques  do  not  always  lead  to  a  design  that  isstableand  atthe 
same  ti  me  robust  to  parameter  uncertai  nti  es. 

1 1  appears  that  a  key  reason  for  i  nstabi  I  ity  at  hypersoni  c  speeds  i  s  that  the  centre  of  pressure 
Cp  position  (without  active  control)  remains  constant  irrespective  of  the  angleof  attack,  M# 
and  altitude.  Hence,  stability  augmentation  due  to  moment  arm  effect  through  shift  in  Cp 
position  (as  in  thecaseof  subsonic  and  low  supersonic  speeds)  does  not  occur  at  hypersonic 
speeds.  Thus  without  activecontrol  a  hypersonic  vehideis  likely  to  belightly  damped  or  even 
unstable.  A  dive  control  provides  for  stability  and  rapid  damping  of  thetransients  (following 
a  disturbance)  by  changing  the  Cp  position  in  a  controlled  fashion.  In  fad,  adivecontrol  can 
often  make  a  vehicle  stable  even  when  the  static  stability  margin  indicates  instability. 

With  uncertain  aerodynamics  and  the  fad  that  Cp  position  remains  fixed,  it  is  clear  that  a 
hypersonic  vehicle  would  requirean  adiveand  robust  (insensitive  to  paramder  uncertainties) 
control  to  accompl  i  sh  stabl  e sustai  ned  fl  i  ght.  M  oreover,  i f  a  hypersoni  c  vehi  cl  e  i  s  expeded  to 
demonstrate  desirable  flight  qualities  over  a  large  flight  envelope  (in  terms  of  M#  and 
altitude)  then  an  adaptive  control  scheme  needs  to  be  implemented. 

Because  of  the  engineering  difficulty  and  performance  costs  associated  with  conventional 
aerodynamic  control  of  hypersonic  vehicles,  several  alternative  approaches  have  been 
proposed  based  on  solely  internal  systems.  One  of  these  is  a  moving  mass  intended  to  effed 
control  by  shifting  the  centre  of  gravity  (Cc)  position  axially  and  laterally.  However,  the 
extent,  speed  and  precision  of  movements  of  the  weights  may  not  be  sufficient  to  confer 
stabi  I  ity  nor  i  ndeed  provi  de  the  control  requi  red  to  mai  ntai  n  the  angl  e  of  attack  at  a  val  ue  to 
achieve  the  desi  red  lift  force.  Moreover,  precision  control  of  angleof  attack  is  necessary  for 
many  high  speed  air-breathing  propulsion  systems  to  operate  satisfadori I y. 

In  this  report  static  stabi  I  ity  and  Cp  position  issues  are  ad  dressed  through  simpleanalysisto 
highlight  some  of  the  key  features  of  hypersonic  aerodynamics.  The  inadequacy  of 
conventional  control  design  techniques  for  a  hypersonic  vehi  decontrol  is  pointed  out.  Also 
included  is  a  sedion  on  control  effedors  (fins,  canards,  thrust  vedoring  etc)  that  are 
commonly  used  in  conventional  missiles  and  have  been  proposed  for  a  number  of  future 
hypersonic  vehicles. 

Analysis  and  design  of  a  hypersonic  vehicle  with  adivecontrol  is  based  on  the  longitudinal 
and/or  lateral  dynamics  model.  A  longitudinal  plane  dynamic  model  for  a  hypersonic 
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airframe  has  been  derived  in  this  report  along  with  alinearised  (small  perturbation)  model  in 
state  space  form.  This  model  can  be  extended  to  indudestructural  vibration  effects  aswell  as 
interaction  with  engine  dynamics  for  an  air  breathing  propulsion  system. 


2.  Aerodynamic  Static  Stability 

Accordingto(Etkin  Band  Reid  LD;  1996),  "An airplanecancontinuein steady  un-accelerated 
fl  i  ght  only  when  the  resultant  external  forces  and  moments  about  CG  both  vani  sh.  Thi  s  i  s  the 
condition  of  longitudinal  balance.  If  the  pitching  moment  were  not  zero,  theairplanewould 
experience  a  rotational  acceleration  component  in  the  direction  of  the  unbalanced  moment." 
This  situation  is  depicted  in  Figure! 


Figurel:  A  erodynamic  forces  &  moments  for  steady  sustained  flight 


Traditionally  'static  margin'  has  been  used  as  a  measure  of  an  air  vehicle's  stability.  This 
criterion  implies  that  the  centre  of  pressure  Cp  lie  behind  (as  measured  from  the  nose)  the 
centre  of  gravity  CG .  Itshould,  however,  be  noted  thatwhilethiswill  givestability  itwill  not 

guarantee  steady  sustained  (or  level)  flight.  The  latter  requires  that  there  should  be  a  positive 
angle  of  attacka,  and  that  this  should  be  maintained  during  a  steady  flight,  in  order  to 
achieve  lift  to  balance  the  weight  of  the  vehicle.  This  condition  known  as  'trim'  typically 
requires  the  application  of  a  control  force  through  an  appropriate  fin  deflection,  creating  a 
positive  angle  of  attack  and  a  shift  in  Cp  position  to  a  location  relative  to  CG  suchthatthe 
resultant  moment  about  CG  is  zero.  If  Cp  is  above  CG  then  it  can  be  shown  that  steady 
sustai  ned  fl  i ght  may  be achi  eved  with  Cp  behi  nd  CG  as  show n  i  n  Fi  gure  2  (compare thi  s  w ith 
Figure  1). 

The  position  of  Cp  (defined  as  the  point  at  which  total  lift  and  drag  forces  act)  does  not 
remain  fixed.  lnfact(Garnell  Pand  East  DJ,  1977)  notes  that:  "It  is  seen  that  at  subsonic  speeds 
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and  very  low  supersonic  speeds  the  extends  to  be  rather  more  forward  than  at  higher  Mach 
numbers.  Also  the  changes  in Cp  with  incidence  can  be  considerable  at  low  speeds;  this  is 
mainly  dueto  general  rearwards  shift  in  the  Cp  of  the  body  with  increasing  incidence,  the 
Cp  of  the  control  surfaces  and  wings  changing  very  little". 


A  plot  of  the  Cp  position  as  a  function  of  the  angle  of  attack  and  M  #,  for  a  fin  controlled 
missile,  isgiven  in  Figure3.TheCp  position,  moves  aft  when  angle  of  attack  a  increases,  and 
moves  forward  when  the  angl  e  of  attack  a  decreases;  i  n  additi  on,  the  I  ift  force  al  so  i  ncreases 
with  increasing  a  .This  produces  a  restoring  moment  that  counters  the  increase  in  a  and  is 
the  very  phenomenon  that  makes  stabl  esustai  ned  fl  i  ght  possi  bl  e.  Thi  s  shift  i  n  the  Cp  positi  on 
as  a  function  of  a ,  occurs  mainly  in  the  subsonic  and  transonic  regimes  but  is  negligi  bleat 
high  supersonic  and  hypersonic  speeds. 


Figure  4  shows  a  typical  plot  for  of  the  pitching  moment  coefficient,  about  Ccfor  a  fixed 

( 1 

elevator,  versus  the  angle  of  attack  « ;  here:  (m  =  CMq*  Sc) ,  and  q*  =  -pv2  is  the  dynamic 

\2 

pressure.  Note  also,  that  the  lift  force  is  given  by:  (L  =  CLq*s) ,  and  the  angle  of  attack  is 
measured  w.r.tto  thezero  lift  line  of  the  vehicle. 
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For  this  case  if  a  disturbance  causes  a  to  suddenly  increase  from  its  nominal  value  a, y  then 
the  negative  moment  (due  to  aft  movement  of  C^and  an  increase  in  CL)  -  a  condition  for 
positive  stiffness  acts  to  restore  a  to  its  nomi  nal  val  ue.  Si  mi  I  arl  y  a  sudden  di  sturbance  i  n  the 
op posi  te  d  i  recti  on  i  nd  uces  a  posi  ti  ve  moment  (d  ueto  forw ard  movement  of  Cp  and  a  decrease 
inCz )  that  acts  to  restore  stability.  It  is  assumed  that  the  velocity  of  the  vehicle  remains 
unaltered  during  these  transient  disturbances.  Key  conditions  for  stability  are: 


-O'  C  =  <  Q 

da 


(i) 


The  moment  i  s  taken  to  be  posi  ti  ve  i  n  the  d  i  recti  on  show  n  i  n  Fi  gu  re  1,  and  a  =  a0  defi  nes  the 
equilibrium  angle  of  attack  for  that  steady  flight  condition.  Note,  also,  that  asymptotically 
stable  flight  is  achieved,  following  a  disturbance,  provided  there  is  sufficient  aerodynamic 
damping. 
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We  now  consider  the  role  that  the  centre  of  pressure  plays  during  thetransient  behaviour  of 
the  vehicle.  In  figure  5,  the  positions  of  Cp  and  CG  are  shown  measured  in  terms  of  their 

respective  distances  xcp,xCGfrom  vehicle's  nose.  The  moment  (about  CG)  equation  may  be 
written  as  (for  small  a ): 

M  =  M  0  —  L(  xcp  —  xCG )  (2) 

Or  equivalently,  in  terms  of  the  moment  and  lift  coefficients  CM,CM0,CL  respectively, 
equation  (2)  may  be  written  as: 

C M  —  CM0  C ,  ( xcp  —  xCG  J  (3) 

c 

CM0  is  the  pitching  moment  independent  of  a.  For  M#  less  than  approximately  5 
(McCormick  BW,  1994;  MieleA,  1962),  CL  =CLa.a  .Thus: 


CM  -  C MO  C  Lu-a\  XCP  ~XCG  j 

c 

For  steady  flight  condition  a  =  a0;  xcp  =  xcpo;  CM  =  0  ,  that  is: 


(4) 

(5) 


xSM  —  ( xc/>  —  XCG )  wil1  be  referred  to  as  the  'static  margin'.  Partial  differentiation  of 
equation  (4)  w.r.t  a ,  we  get: 
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a.  Thus,  if  — —  =  0  (e.g.  hypersonic  flight),  then  conditionsfor  stability  imply  that: 

da 

xcp  >  xcg  >  -e.  CP  must  be  aft  of  CG  for  the  vehi  cl  e  to  remai  n  stabl  e. 

dx 

b.  How  ever,  i  f  — —  >  0 ,  (e.  g.  su  bson  i  cor  transon  i  c  f  I  i  ght  or  acti  ve  control  I  ed  f  I  i  ght), 

da 

dx 

then  conditionsfor  stability  imply  that  xcp  >xCG  -  a — —  .  I  nth  is  case  steady  and 

da 

stabl  e  f I  i  ght  i  s  ach  i  evabl  e  w  i  th  zero  or  negati  ve  stati  c  margi  n . 
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3.  Effect  on  Stability  Due  to  Shift  in  CG  Position 

It  was  noted  in  Section  2  that  as  the  M#  increases  the  shift  in  the  Cp  due  to  changes  in  the 

angle  of  attack  becomes  progressively  smaller  and  isnegligibleathigh  M#.  (Ostapenko  NA, 
1980)  notes  that  for  high  supersonic  flow:  "The  theoretical  investigation  of  the  aerodynamic 
characteristics  of  circular  cones  has  shown  that  their  centre  of  pressure  does  not  depend  on 
the  angl  e  of  attack  w  hen  the  shock  wave  i  s  attached  to  the  apex  of  the  cone.  1 1  was  establ  i  shed 
experimentally  for  star-shaped  bodies  that  the  position  of  the  centre  of  pressure  for  such 
bod  i  es  hard  I  y  changes  i  n  a  w  i  de  range  of  M  ach  numbers  and  angl  es  of  attack" .  I  n  vi  ew  of  thi  s 
and  other  published  research,  it  will  be  assumed  that  Cp  position  remains  constant  at 

hypersoni  c  speeds  irrespective  of  M#and  angle  of  incidence. 

Let  us  now  consider  the  case  where  the  exposition  suddenly  moves  aft,  causing  a  positive 
moment  to  occu  r  and  the  vehi  cl  e  angl  e  of  attack  to  i  ncrease  from  i  ts  equ  i  I  i  bri  um  val  uea0 .  A t 
low  M#  the  aft  movement  of  C^and  increasein  CL  acts  to  restore  vehicle  stability  about  a 
new  equilibrium  angle  of  attack  a,  .However,  at  high  M#theaft  shift  of  Cp  position  does  not 
occur  and  increase  in  C,  may  not  be  sufficient  to  restore  stability.  In  this  case  the  vehicle 
cou  I  d  become  u nstabl  e  and  toppl  e over  u n I  ess  acti  ve  control  s  are avai  I  abl  e.  A  cti  ve control  al  so 
allows  an  air  vehicle  to  achieve  new  trim  conditions  to  counter  changes  in  CG. 

Oneobvious  way  to  avoid  instability  is  to  design  thevehicle  with  CG  significantly  forward  of 
the  Cp  (large  static  margin);  however  thiswould  makethevehidesluggish  in  responding  to 
demands  in  attitude  and  flight  path  changes  -  a  characteristic  not  very  desirablefor  some 
applications.  Even  with  large  static  margin,  however,  activecontrol  will  still  be  required  to 
achieve  steady  sustained  flight. 

U  n  I  i  ke  su  bson  i  c  and  su  personi  c  veh  i  cl  es,  the  stabi  I  ity  of  hy  person  i  c  veh  i  cl  es  ( uti  I  i  si  ng  acti  ve 
control)  cannot  be  assessed  purely  from  static  margin  alone.  In  fact,  as  stated  in  a  study  by 
(Johnson  DB,  Thomas  R,  and  Manor  D;  2001):  "Static  margin  has  been  the  standard  indicator 
of  I  ongi  tud  i  nal  stabi  I  ity  for  many  years.  H  owever,  experi  ence  w  ith  previ  ous  ai  r  vehi  cl  es  fl  own 
with  si  gnifi  cant  I  evel  s  of  stati  c  i  nstabi  I  ity  i  ndi  cated  that  the  con  venti  onal  stati  c  margi  n  i  s  not  a 
valid  indicator  for  these  vehicles.  In  one  case,  30%  unstable  was  flyable,  whilein  another  15% 
unstable  was  completely  unacceptable".  That  is,  in  certain  cases,  the  hypersonic  vehicle  with 
activecontrol  wasstablealthoughthestatiemargin  indicated  instability. 
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4.  Hypersonic  Flight  Stability  and  Control  Issues 


Hypersonic  vehicle  aerodynamics  and  its  longitudinal  stability  and  control  characteristics 
have  been  studied  by  a  number  of  authors.  For  example:  aerodynamic  flow  characteristics  of  a 
hypersoni  c  glidevehi  cl  e(M#>5)  was  considered  by '(MieleA;  1962)'.  The  main  characteristics 
of  this  flow  are  that: 

a.  The  shock  waves  originating  attheleading  edge  of  the  body  lie  close  to  the  body  so 
that  the  interaction  with  the  body  is  strong. 

b.  Hightemperaturesexistintheregionsbetweentheshockwavesandthebody  and  it 
may  be  necessary  to  consider  real  gas  effects  (molecular  vibration,  dissociation,  and 
ionisation)  when  analysing  theflow  fields. 

c.  At  very  hi  gh  M  #,  theshock  waves  may  be  assumed  to  be  al  most  identi  cal  to  the  body, 
at  I  east  at  the  front  portion  of  the  body,  and  the  molecules  crossing  theshock  waves 
conserve  the  tangential  component  of  the  velocity  but  lose  most  of  the  normal 
component. 


In  view  of  the  above,  a  possible  design  of  a  glide  vehicle  operating  at  high  M  #would  haveto 
be  a  compromise  between  aerodynamic  and  heat  transfer  requirements.  For  example,  the 
glide  vehicle  could  have  lifting  surfaces  with  planar  edges  and  a  nose  that  is  blunt. 
Aerodynamic  lift  and  drag  coefficients:  CL,CD ;  zero  lift  drag  coefficient  CD0 ;  aerodynamic 


efficiency  Emax 


max 


;  the  induced  drag  factor  K ,  and  the  exponent  n  (as  defined  in 


equation  10 below)  aredepicted  in  Figure 6. 


It  is  clear  that  key  differences  exist  between  subsonic/  supersonic  and  hypersonic 
aerodynamic  characteristics.  For  example,  for  a  supersonic  vehicle  1.2<  M  #<5,  the  lift  and 
drag  coefficients  are  given  by  '(Etkin  B  and  Reid  LD;  1996)': 


C,=CLaa  (9) 

CD  -  cDo  +  KCl"  (10) 


The  val  ue  of  n  =  2  for  M  #<5,  and  the  three  constants 

Cl,CD0,K  are  functions  of  the  configuration,  thrust  coefficient  and  M#.  Whereas  for  a 


3 

hypersonic  vehicle  M#>5,  «  =  — ,weget: 

( l  .  A 

CL=  —  CNa  sin  2a  +  CNaa  sin  a  sin  a  cos  a 

\2  ) 

3 

(11) 

CD=CDo+KCL2 

(12) 

8 


DSTO-TR-2358 


Here  CNa  =(cLa\a=0),CNatt  are  coefficients  (independent  of  a)  dependent  on  M#and 
configuration.  For  0  <  a  <  180 ° : 


CLa  =  CNa  (cos1  a-  2  sin2  a  cos  a ) 

+  CNaa  ( 2  sin  a  cos 2  a- sin3  a} 


(13) 


For  small  values  of  the  angle  of  attack  a ,  equations  (11)-(13)  may  be  written  as: 

CL=(CNaa  +  CN<ma2)  (14) 

CLa=CNa(l-2a2)  +  CNaa(2a-a3)  (15) 


Figure6:  A  erodynamic  characteristics  of  a  hypervelocity  glider  ('M  iele;  1962') 


It  is  also  noteworthy  (MieleA;  1962)  that  the  maximum  value  of  the  ratio  L/D  isequal  to 
5-10  for  supersonic  vehicles  and  only  about  1-5  for  hypersonic  vehicles.  Also,  unlike 
supersoni  c  vehi cl  es  the  I  ift  and  drag  coeffi  ci ent  remai  n  constant  for  hypersoni  c  vehi cl  es  for 
M  #>6. 

A  numerical  example  (Etkin  B,  1972),  with  parameters  given  in  Appendix  1,  is  used  to 
demonstrate  typical  forces  and  moments  characteristics  of  a  hypersonic  vehicle.  Values  of 
CL,CD,CL  ,CD  for  this  example  are  plotted  in  Figures  7(a)-(d)  for  both  exact  and  approximate 

val  ues  of  the  vari  ous  parameters;  i  t  i  s  seen  that  the  approxi  mati  on  i  s  val  i  d  for  smal  I  val  ues  of 
angl  e  of  attack  a  <12° . 
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Static  stability  analysis  can  be  conducted  for  this  example;  since  CL  is  now  a  nonlinear 
function  of  a ,  equations  (3)  -  (7)  may  be  used  to  compute  the  static  stiffness  curve.  Ignoring 
the  shift  in  the  (^(the  last  term  in  equation  (7)),  the  following  values  are  obtained  (see 
Appendix  1): 


a,  (deg) 

XSM 

c mo  (N  m) 

CMa  (N  m) 

10.0 

1.445 

0.00493 

-0.095  CLa(a) 

Fi  gure  8  i  s  the  pi  ot  of  stati  c  stabi  I  ity  characteri  sti  c  of  the  vehi  cl  e  i  n  the  exampl  e;  it  shows  that 
the  vehicle  is  statically  stable  with  positive  stiffness. 

While  the  above  analysis  might  suggest  that  the  vehicle  will  be  stable  without  any  active 
controls,  however,  aerodynamic  design  based  on  static  analysis  aloneisunreliableand  a  more 
comprehensive  analysis  is  necessary  using  the  dynamic  model  for  hypersonic  vehicle 
considered  later  in  this  report. 

Several  authors  have  reported  thatthereis  significant  disparity  bdweenthewind  tunnel  data 
(aerodynamic  coefficients)  and  predictions  based  on  this,  as  compared  totheactual  flight  data 
(I  I  iff  KW  and  Shafer  MF;  1993).  Figure  9  shows  a  plot  of  the  actual  (obtained  from  flight 
tests)  exposition  and  the  predicted  values  (obtained  from  ground  tests).  Clearly  in  order  to 
ensure  acceptable  aerodynamic  performance  of  a  hypersonic  vehicle,  there  is  a  need  to  have 
an  activecontrol  and  this  has  to  be  robust  enough  to  cater  for  uncertainties  in  the  predicted 
and/  or  computed  aerodynamic  design  parameters. 

A  control  system  based  on  moving  the  centre  of  gravity  (i.e.  moving  internal  weights)  still 
uses  lift  and  drag  forces  to  produce  the  required  control  moments.  Flowever,  it  may  be 
diffi  cult  to  achi  eve  rapi  d  changes  i  n  the  moments  that  are  needed  to  gi  vefast  enough  control 
response.  The  use  of  this  mechanism  to  achieve  stable  and  steady  sustained  flight  at 
hypersonic  speeds,  apart  from  obvious  engineering  difficulties  may  not  necessarily  achieve 
the  demonstrati  on  objective.  In  fact,  various  reported  studies  on  stability  and  control  indicate 
that  a  hypersonic  vehicle  tend  to  beunstableand  require  active  control  to  achi  eve  desirable 
transient  properties  (McLaen  D  et.al.;  2007). 
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F  i gu  re  7a:  E xact  an  d  approx i mate  ( seri es  2)  I i  ft  coeffi ci en  t  CL 


Exact/Approx.  CD 

.0 

i  » 

- - - - 

0  - 

0 

0  2 

0  4.0  6 

0  8.0  1C 

Angle  of  Attack 

Seriesl  — ♦— Series2 

.0  12 

.0  14 

.0  16 

Figure  7b:  Exact  and  approximate  (series  2)  drag  coefficient  CD 
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F  i gu  re  7c:  E xact  an  d  approximate  ( seri es  2)  the  coeffi ci en  t  CLa 


Figure  Id:  Exact  and  approximate  (series  2)  the  coefficient  CDa 
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Figure8:  Static  stability  plots  CM  against  a 


Figure9:  TheSTS  longitudinal  aerodynamic  Centre  of  P  ressure  location  comparison  (‘RomerePO 
and  M  ilesAW,  1983’) 
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5.  Hypersonic  Aerodynamic  Model  forControl  Analysis 

and  Synthesis 

I  n  this  section  weconsider  the  longitudinal  dynamical  model  for  a  hypersonic  vehidethat  is 
employed  for  stability  and  control  analysis  and  synthesis.  A  small  perturbation  model 
(linearised  model)  is  also  considered.  These  models  may  be  used  to  study  the  stability 
behaviour  for  changing  flight  conditions  and  for  control  sensitivity  and  robustness  studies. 


The  longitudinal  force  and  moment  equations  (see  Figure  10)  for  a  hypersonic  vehicle  are 
derived  in  Appendix  -1,  equations  (A-2.10)-(A-2.14),  these  are  reproduced  here,  for 
convenience.  These  equations  define  the  non-linear  longitudinal  model  that  is  used  for 
analysis  and  synthesis  of  stability  and  control  of  a  hypersonic  vehicle: 


— V - gsin(O-a) 

dt  m  V  ’ 

d  Tzw~L  gcos(d -a) 

— a- — 5 - b - - - L-io  +  q 

dt  mV  V 

d  M 
— 9-  — 

*  K 

d  e  V  cos(0-a) 

— 6  =q  +  io  + - - - - 

dt  R 


—  R  =  V  sin(d-a) 
dt  V  ’ 


(16) 

(17) 

(18) 

(19) 

(20) 


14 


DSTO-TR-2358 


These  equations  are  similar  to  those  presented  by  (Bolander  M,  2009);  (Bilimora  K  and 
Schmidt  D;  1995;  Groves  KP,  et.  al 2005)  except  for  the  definition  of  the  pitching  moment  q 
used  inthecurrent report; thepitching moment  Q  used  bytheaboveauthorsisequivalentto: 
Q  =  q  +  (oe  +//.  These  authors  neglect  earth's  rotation  and  also  assume  that//  =  0  . 

The  linearised  longitudinal  model  (small  perturbation  model)  has  also  been  derived  in 
Appendix- 1,  equation  (A-31),  and  may  beused  for  steady-state  stability  analysisand  control 
system  synthesis.  Both  the  non-linear  and  the  linearised  model  require  knowledge  of  the 
aerodynamic  parameters  (lift,  drag,  and  moment  coefficients).  These  parameters  can  be 
obtai  ned  either  vi  a  theoreti  cal  model  s  (A  nderson  J  D ,  2006)  or  experi  mental  I  y  through  grou nd 
(wi  nd  tunnel )  testi  ng.  Because  of  the  I  ack  of  avai  I  abi  I  ity  of  actual  fl  ight  data,  the  parameters 
used  inthemodel  are  predicted  values  and  may  havesignificantuncertainty  attached  to  them. 
Hence  the  various  techniques  that  have  been  suggested  for  control  system  design  generally 
uti  I  i  se  robust  or  adapti  ve  control  techni  ques  to  al  I  ow  for  parameter  vari  ati  ons  and  uncertai  nty 
(Faruqi  FA  and  Jijoong  K,  2009; Jan kovsky  P  et.  al.  2007). 

H  ypersoni  c  fl  i  ght  control  requi  res  a  hi  gh  degree  of  preci  si  on,  and  i  n  many  cases  even  greater 
than  for  a  conventional  aircraft.  In  addition,  the  control  system  needs  to  allow  for  physical 
effects  such  as  heati  ng  that  may  cause  damage  to  wi  ri  ng  and  el  ectroni  cs,  sensor  degrad  ati  on, 
structural  distortion  and  erosion  of  control  surfaces.  Methodsof  cooling  and  insulation  have 
been  suggested  that  may  mitigate  these  problems.  I  n  the  case  of  an  air-breathing  (Ramjet  or 
Scram-jet)  moment  and  force  interaction  between  the  engine  and  airframe  becomes  an 
important  consideration  and  methods  have  been  proposed  for  including  these  effects  in  the 
vehicle  dynamics  (Bolander  M,  2009). 

From  a  navigation  and  guidance  perspective,  our  past  experience  suggests  that  for  a 
hypersonic  vehicle  (as  for  the  case  of  a  supersonic  or  subsonic  vehicle),  a  suite  of  inertial 
sensors  (I MU),  GPS  and  other  sensors  (e.g.  magnetometers,  visual  imaging  sensors)  will  be 
required  to  obtain  data,  such  as:  position,  velocity  and  attitudefor  navigation.  Ifthevehideis 
required  to  guide  to  a  target  then  seeker  data,  such  as:  target  rangeand  range  rate,  target  LOS 
angleand  rates,  lookangleand  rates  areal  so  required.  Issues  that  must  be  addressed  include: 
sensor  accuracy,  data  rate,  sensor  placement  in  theairframe  and  airframe  structural  vibration 
and  heating  effects.  Thelatter  aspect  is  particularly  relevant  to  seeker  radomedesign. 


6.  Control  Effectors  Options 

Hypersonic  airframe  configuration  and  design  requires  an  understanding  of  the  flow-fields 
that  exist  around  the  vehicle,  the  heat  transfer,  and  stability  and  control  issues.  In  addition, 
si  nee  hi  gh  I  i ft/  d rag  rati  os  are  not  achi  evabl  e  at  hypersoni  c  speeds,  a  parameter  opti  mi  sati  on 
approach  is  generally  used  to  identify  the  desired  airframe  and  control  configuration 
(Bowcutt,  KG.,  Anderson,  J  D.,  and  Capriotti,  D,  1987).  Most  of  the  reported  designs  of 
hypersonic  vehicles  have  been  confined  to  the  wave-rider  and  this  has  incorporated  rudder 
and  combined  elevator  and  aileron  (at  the  rear)  as  control  surfaces  (Cockrell  CE  et.al.,  1995). 
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When  designing  a  hypersonic  vehicle,  however,  it  is  necessary  to  consider  various  other 
control  effector  options  such  as  canards  and/  or  thrust  vectors.  For  example:  proposed  control 
method  for  H yStri ke  (H  igh  Speed  Strike  Missile)  and  Fast  Flawk  (M#4,  Low  Cost  Missile 
System  -  LCMS)  isTVC  achieved  through  body  bending  (Figure  11);  the  FI  i Fly  (Flypersonic 
Flight  Demonstrator)  and  Falcon  (FI  ypersonicCruiseVehide)  use  tail  control  (Figures  12  and 
14);  theShyfe  (Sustained  Flypersonic  Flight  Experiment)  uses  canard  for  longitudinal  and 
lateral  control  (Figure  13).  Boeing's  proposed  X51  (Flypersonic  Weapons  Technology 
Demonstrator)  appears  to  use  tail  control  (Figure.  15). 

Tail  control  is  probably  the  most  commonly  used  form  of  conventional  (subsonic  and 
supersonic)  missile  control,  particularly  for  longer  range  applications,  and  appears  in  a 
number  of  hypersonic  mi  ssi  I  econcepts.  This  is  mainly  because  tail  control  provides  excel  lent 
manoeuvrability  at  the  high  angles  of  attack,  which  are  needed  to  intercept  manoeuvrable 
aircraft.  Missiles  using  tail  control  may  also  befitted  with  fixed  wings  or  strakes  to  increase 
lift  and  improve  range. 

Canard  control  has  been  used  in  some  conventional  missiles,  particularly  short-range  air-to- 
ai  r  mi  ssi  I  es,  and  has  been  proposed  for  Shyfe.  The  key  advantage  of  canard  control  is  that  it 
provides  better  manoeuvrability  at  low  angles  of  attack,  but  canards  tend  to  become 
i  neffecti  ve  at  hi  gh  angl  es  of  attack  because  of  fl  ow  separati  on  that  causes  the  canard  to  stal  I . 
Si  nee  canards  are  ahead  of  the  centre  of  gravity,  they  cause  a  destabi  I  i  si  ng  effect  and  requi  re 
large  fixed  tails  to  keep  the  missile  stable.  The  fins  usually  provide  sufficient  lift  to  make 
wings  unnecessary.  One  way  of  avoiding  canard  from  stalling  is  to  have  double  or  split 
canard  (two  sets  of  canards  in  close  proximity,  usually  one  immediately  behind  the  other).  In 
the  split  canard  configurati  on, thefirst  canard  setisfixed  whilethesecond  set  is  movable.  The 
first  set  of  canards  generates  strong,  vortices  that  increase  the  speed  of  the  airflow  over  the 
second  set  of  canards  making  them  more  effective. 


21”  DIA  X  256  LONG  -  3152  LBM  (SURFACE  LAUNCHED) 
21”  DIA  X  168”  LONG  -  1880  LBM  (AIR  LAUNCHED) 


Figurell:  HyStrike  (H  igh  Speed  Strike M  issile) 
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Figure  12:  FiiFly  (FI ypersonic  Flight) 
D  em  on  st  rat  or) 


I  n  ad  d  i  ti  on ,  the  vorti  ces  d  el  ay  f  I  ow  separati  on  and  al  I  ow  the  canard  s  to  reach  h  i  gher  angl  es  of 
attack  before  stalling.  The  canards  also  tend  to  interact  with  the  fins  and  introduce  induced 
roll  which  requires  additional  roll  control  installed  in  the  fins. 

W  i  ng  control  has  been  used  i  n  some  of  the  earl  i  er  conventi  onal  mi  ssi  I  e  but  i  t  i  s  becomi  ng  I  ess 
common.  Most  missiles  using  wing  control  are  longer-range.  The  key  advantage  of  wing 
control  isthatthedeflectionsofthewingsproducea  very  fast  responsewith  littlemotion  of 
the  body.  Thi  s  feature  results  i  n  smal  I  seeker  tracki  ng  error  and  al  I  ows  the  mi  ssi  I  e  to  remai  n 
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locked  on  target  during  large  manoeuvres.  The  mai n  disadvantage  is  that  the  wings  must 
usual  I  y  be  quite  I  arge  i  n  order  to  generate  both  suffi  ci  ent  I  ift  and  control  effectiveness,  w  hi  ch 
makes  the  missiles  rather  large  overall.  In  addition,  the  wings  generate  strong  vortices  that 
i  nteract  w i th  the  f i  ns  cau si  ng  the  mi  ssi  I  e  to  rol  I . 

Th  ru  st  vectori  ng  i  s  a  method  of  d  ef  I  ecti  ng  the  mi  ssi  I  e  exhau  st  to  generate  forces  and  moment 
to  cause  the  vehicle  to  manoeuvre  in  longitudinal  and  lateral  direction.  Another  technique 
that  is  currently  being  investigated  is  the  so  called  reaction  jets  control.  Reaction  jets  are 
usually  small  ports  in  the  surface  of  a  missile  that  create  a  jet  exhaust  perpendicular  to  the 
vehi  cl  e  surface  and  produce  an  effect  si  mi  I  ar  to  thrust  vectori  ng.  These  techni  ques  give  hi  gh 
off-bore-sight  capability  to  conventional  missiles.  The  key  advantage  of  such  controls  is  that 
they  can  function  at  very  low  speeds  or  in  a  vacuum  where  there  is  littleor  no  airflow  to  act 
on  conventional  fins.  Theprimary  drawback,  however,  isthat  they  will  not  fundi  on  once  the 
fuel  supply  isexhausted.  Note  that  most  missiles  equipped  with  readion  jets  controls  do  not 
rely  on  these  controls  alone  for  manoeuvrability,  but  only  as  a  supplement  to  aerodynamic 
surfaces  I  i ke  canards  and  tai I  fi  ns. 

Finally,  it  should  be  noted  that  some  missiles  have  in  the  past  used  conventional  controls 
si  mi  lar  to  those  empl  oyed  by  ai  rcraft.  These  systems  are  usual  I  y  referred  to  as  bank-to-turn 
control  s  si  nee  the  mi  ssi  I  e  banks  much  I  i  ke  an  ai  rpl  ane  wou  I  d .  A  nother  form  of  control  s  that  i  s 
currently  under  investigation  is  the  use  of  micro-flow  that  alters  the  flow  field  around  the 
ai  rframe  of  the  vehicle  and  generates  forces  and  moments  sufficient  to  turn  the  vehicle  i  n  a 
controlled  fashion.  Insufficient  data  is  currently  available  on  this  to  merit  serious 
consideration.1 


7.  Conclusions 

This  report  has  been  produced  in  order  to  address  aerodynamic  charaderi sties  and  other 
i  ssues  rel  ati  ng  to  hy  personi  c  vehi  cl  es  that  are  deemed  to  be  si  gnifi  cantl  y  d  iff  er  ent  from  those 
of  the  conventional  (subsonic  or  supersonic)  air  vehicles,  guided  missiles  and  air-borne 
weapons.  In  particular  we  have  addressed  issues  that  are  relevant  to  stability  and  control  of 
hypersonic  vehicles.  This  report  should  add  to  the  existing  knowledge  and  experience  of 
mi  ssi  I  e  gui  dance  and  control  engi  neers  and  make,  other  researchers  and  engi  neers  i  nvol  ved  i  n 
hypersonic  experimentation,  aware  that  these  vehicles  may  not  be  dynamically  stable  and 
require  adive  control  augmentation  in  order  to  achieve  and  maintain  desirable  flight 
charaderi  sties. 


1  Acknowledgement:  Sedions  35  -  39  have  relied  heavily  on  the  article  by  Jeff  Scott  at 
Aerospaceweb.org  Ask  Us- M  issile Control  Systems 
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From  the  reported  work  on  hypersonicvehidestability  and  control  and  some  simple  analysis 
presented  i  n  thi  s  report,  we  are  abl  eto  hi  ghl  i  ght  aerody  nami  c  characteri  sti  cs  that  may  have  a 
strong  bearing  on  our  approach  to  the  analysis,  synthesis  and  performance  evaluation  of 
hypersoni c  weapons  systems.  The  key  fi  ndi  ngs  and  recommendati  ons  are  as  fol  I  ows: 

a.  Unlike  subsonic  and  low  supersonic  vehicles,  centre  of  pressure  position  does  not 
change  for  hypersonic  vehicles  for  changes  in  M#,  angleof  attack  and  altitude.  Active 
control  will  be  required  to  maintain  stability  particularly  in  case  of  changes  to  the  CG 
position  (duetofuel  burn-off)  and  maintain  adequate  damping  during  flight. 

b.  It  has  been  noted  that  aerodynamic  parameters  derived  from  ground  tests  and 
theoreti  cal  consi  derati  ons  do  not  refl  ect  the  actual  i  n-fl  i  ght  aerodynami  cs  at  hypersonic 
speeds.  Uncertainties  in  these  parameters  require  a  control  system  design  based  on 
robust  techniques.  Also,  if  the  flight  envelope  (M#,  altitude)  is  expected  to  vary  then 
adaptive  techniques  are  required  to  maintain  desirable  vehicle  flight  performance. 

c.  At  hypersonic  speeds  thedrag  and  I  iftforces  become  non-linear  functions  of  theangle 
of  attack.  At  the  same  time,  following  a  transition  phase  during  the  transonic  and 
supersonic  speeds,  the  drag  and  lift  coefficients  attain  constant  values  at  hypersonic 
speeds. 

d.  Compared  to  subsonic  and  supersonic  speeds,  the  maxi  mum  value  of  lift/  drag  ratio 
for  hypersonic  vehicles  is  significantly  lower.  For  example,  the  maxi  mum  value  of  the 
ratio  L/D  is  equal  to  5-10 for  supersonic  vehicles  and  only  about  1-5  for  hypersonic 
vehicles.  Research  is  on-going  to  exploreways  of  increasing  this  ratio. 

e.  Heating  effects  are  significant  at  hypersonic  speeds  and  it  has  been  suggested  that 
rounding  of  the  nose  and  other  leading  edges  may  be  required  to  reduce  thermal 
gradients.  Ablative  material  to  dissipate  heat  has  also  been  suggested.  H  eating  has  a 
negative  affect  on  vehicle  structural  integrity  and,  depending  on  the  airframe,  may 
cause  structural  vibration  which  needsto  be  catered  for  by  appropriate  control  system 
design. 

A  longitudinal  plane  dynamic  model  aswell  asasmall  perturbation  model  has  been  derived 
in  this  report.  These  may  beused  for  stability  analysisand  control  system  design.  The  report 
al  so  i  nd  udes  a  summary  of  the  different  types  of  control  effectors  that  have  been  used  i  n  the 
past  on  a  variety  of  supersonic  and  hypersonic  air  vehicles.  Also  included  are  some  well 
known  airframe  configurations  that  have  been  proposed  for  a  number  of  current  hypersonic 
vehicle  programs.  The  object  is  to  present  control  effector  options  that  are  available  to  a 
hypersonic  vehicle  designer.  The  preferred  selection  clearly  depends  on  the  operational 
requirement  and  numerous  engineering  considerations,  and  will  be  dictated  by  a  multi¬ 
disciplinary  team  responsible  for  total  system  design  and  performance  evaluation. 
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Appendix  A:  Hypersonic  Vehicle  Example: 


A  hypothetical  hypersonic  vehicle  is  considered  at  a  nominal  altitude  of  53  km,  andaspeed 
relative  to  earth  of  approx  7.4  km/  s 

Aerodynamic  Parameters 
CL  =1.75  sina.cosa.|sina| 

CD  =0.0625+1.69  sin2a.|sina| 

CLa  =3.5|sina|cos2a  -1.75sin2a|sina| 

CDa  =5.07sina.cosa|sina| 

Cmq  =  -0.03 
Cma  =  -0.055 
a  =  10  degrees 
c  =  15.25m 


Approximate  Values 
CL  =1.75  a2 
CD  =0.0625+1.69  a3 
CLa  =3. 5a  -1.75a3 
Coa  =5. 07a2 

The  val  ues  given  i  n  the  above  tabl  es  are  pi  otted  i  n  Fi  gures  6(a)-(d)  for  both  the  exact  and  the 
approxi  mate  val  ues. 


We  use  equation  (15)  for  calculating  the  static  margi  n  by  setti  ng  a  =  10°  and 


=^-  =  0,  that  is: 
da 


( [XCP0  XCG ) 


c.C. 


CLa(a0 ) 


15.25  x  0.055 
0.5803 


=  1.445m 


CMa=--CLa(a)(xCp0-xCG)^-0.095CLa(a) 


(A- 1.1) 
(A -1.2) 


From  equation  (3),  we  get: 

C M0  —  —C ' 0  )iXd'0  —  XCG  ) 


=  0.00493 


0.05197x1.445 

15.25 


(A- 1.3) 
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Appendix  B:  Hypersonic  Longitudinal  Dynamic  Model 


FromFigure9,theforceand  moment  equationsforthelongitudinal  motion  of  avehideinthe 
flight  path  (wind  axes)  (Etkin  B,  1972)  may  be  written  as: 

Tx  -  D  -  mg  sin  y  -  mV 
Tzw~L  +  ms  cos  y  =  -mV  (qw  +  qeH, ) 

Where: 

Tx  :  is  the  thrust  x-component  of  the  thrust  vector. 

T_  :  i  s  the  thrust  z-component  of  the  thrust  vector. 

D:  is  the  body  drag  force. 

L:  i s  the  body  lift  force. 
m:  is  the  vehicle  mass. 
g:  is  the  earth's  gravity. 

V:  is  the  vehicle  flight  path  velocity. 
qn, :  i  s  the  fl  i  ght  path  rotati  on  rate  about  the  yK  -axes. 

qew:  is  the  earth's  rotation  projected  along  the  yw  -axes. 
q:  is  the  vehicle  pitch  rate  about  the  body  y  -axes. 

M:  is  the  vehicle  pitch  moment  about  the  body  axes. 

/  :  i  s  the  vehicle  moment  of  inertia  about  its  v -axes. 

yy  J 

(Ok,xw,zw):  are  the  x,  z  coordinates  of  the  wind-axes  system  (or  the  flight  path  axes  system) 
with  its  xw  pointing  along  the  vector  V  . 

(Ov  ,xv  ,zv ) :  arethex,y,  z  coordinates  vehicle-axes  system  (defined  astheaxes  system  attached 
tothevehideCG.and  moving  with  it,  with  zv  axis  point  along  thegravity  vector.  The  xv  -axis 
is  taken  to  poi  nt  north.  N  ote  that: 

qw-q-a  ;and0  =  y  +  «  (A-2.4) 


(A-2.1) 

(A-2.2) 

(A-2.3) 


For  the  purpose  of  stability  analysis,  we  may  assume  that  the  vehicle  is  moving  eastward 
along  fixed  latitude^ ,  say  the  equator  (for  the  equator^  =  0 ).  Earth's  rotation  coe  may  be 
included  in  the  vehicle  dynamical  equations  as  follows: 


cos  X 

~1 

e 

= 

0 

ioe  - 

0 

-  sin  X 

0 

(A-2.5) 
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Where: 

cov :  is  the  earth's  rotati  on  vector  w .  r.t  the  vehi  cl  e-axes 

coew  =  [ Pi  9l  rn  J  istheearth's  rotation  vector  w.r.tthewind  axes.  In  thewind-axes system, 
for  the  vehicle  travel  ling  east,  equation  (A-5): 


e  F  rri  W~ 1  6 


0  cosy  -siny 

1 

"  0  " 

-10  0 

0 

coe  - 

-of 

0  siny  cosy 

0 

0 

e  e 

< / ,,  = 

If  we  consider  the  motion  of  the  vehicle  w. r.t  to  earth's  longitude;/ ,  then  weget: 

7=9K-9l+/i^q-a  +  a)e+ju  (A-2.7) 

9  =  y  +a  =  q  +  coe+  ;/  (A-2.8) 


Where: 

=  cosy  R  js  the  location  ofthe  vehicles  mass  centrew.  r.t  earth  and: 

R 

R-Vsiny  (A-2.9) 

Equations (A-2.1)-(A-2.3),  (A-2.7)  and  (A-2.8)  definethevehidedynamics and  may  bewritten 
as: 


1  T  -D 

“  i/  _ 

-  g  sin  ( 0  -  a ) 

(A-2.10) 

V  — 

dt  m 

d  Tzw  L  gcos(9-a ) 

— a- — = - + - - - co  +a 

(A-2.11) 

dt  m  V 

V 

d  M 

— 9  = - 

*  K 

(A-2.12) 

d 

V  cos(O-a) 

(A-2.13) 

— 0  -q  +  co  + 

dt 

R 

—R  =  V  sin(0 
dt  v 

-a) 

(A-2.14) 

These  equations  are  similar  to  those  suggested  by  (Boland  er  M,  2009;  BilimoraK  and  Schmidt 
D;  1995;  Groves  KP,  et.  al„  2005)  except  for  the  definition  ofthepitching  moment/;  used  in  the 
current  report;  the  pitching  moment  Q  used  by  the  above  authors  is  equivalent  to: 
Q-q  +  coe  +  fi .  These  authors  neglect  earth's  rotation  assume  that//  =  0  . 
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Equations  (A-2.10)  -  (A-2.14)  are  non-linear  longitudinal  model  for  a  (hypersonic)  vehicle 
dynamics;  the  system  states  are(V,a,q,  y,R) .  For  small  perturbation  the  above  equations  may 
be  linearised  about  a  given  steady  state;  from  the  above  equations  we  get  the  following 
rel ati onshi  p  for  the  steady  state: 


0  =  ZL_^;(j  T  y  =  o) 

mx  ’ 

0  -  -L0  +  mg0  -  m  V0  (o/  -  q0 ) 

T_  -  -T03a 
0  = 


0=qt+a>‘+jj- 


Equations  (A-2.16)  and  (A-2.19) 

L0=mg0-m0 

The  linearised  longitudinal  model  for  the  vehicle  is  given  by: 


(  y2^ 

2V„(oe+-± 

V. 


(A-2.15) 

(A-2.16) 

(A-2.17) 

(A-2.18) 

(A-2.19) 


(A-2.20) 


d  3T-3D  (  ^ 

—3V  = - g0(30-3a) 

at  m 

—3a  = 
lit 


2of  1 
- 1 - 

R 


V  V, 


o  J 


SV-M «L  +  f*+i, 

mV,  mV,  V, 


il  c  3M 
— 3q  = - 

dt  K 

d  xn  x  dV  Vo  xt, 

— 36  =  3q  + - \3V 

dt  *  ll„  Rj 


~3R  =  V0(36-3a) 


(A-2.21) 

(A-2.22) 

(A-2.23) 

(A-2.24) 

(A-2.25) 


The  aerodynamic  moments  and  forces  (for  hypersonic  flight)  may  bewritten  in  terms  of  the 


respective  coefficients  as: 

3T  =  Tv3V  +  Tr3R  (A-2.26) 

3D  -Da3a  +  Dr3R  +  Dq3q  (A-2.27) 

3L  -La3a  +  Lr3R  +  Lq3q  (A -2. 28) 

3M  =Ma3a+Mq3q  +  Mn3q  (A-2.29) 

3g  -  —3R  -  -2—3R  (A-2.30) 

clR  R0 


Where:  // :  is  the  control  deflection. 
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Combining  equations  (A-2.21)-(A-2.29),  we  may  write  the  linearised  (or  small  perturbation) 
model  as: 


d_ 

dt 


~3V~ 

3a 

3q 

- 

39 

3R 

f  t 

Lv 

f  D  ^ 

a 

So 

[m  J 

l  m  ) 

f  2ioe  1  1 

| 

(  T  +  T  \ 
10  +  « 

V  Vo  Ro  J 

l  mVo  J 

(  \ 

Ma 

(Mo) 

I 

I 

v  yy  / 

V  yy  / 

So 

0 


T  -D 

R  R 

111 


2g0 


mV0  R0V0  j 


fMR"' 


\  J 


(vA 

0 

0 

0 

0 

\R0  J 

VRo) 

0 

-Vo 

0 

Vo 

0 

r  3v' 

r  d  i 

n 

3a 

K 

3q 

+ 

Mo 

39 

0 

[3R 

0 

3,1  (A-23 ]) 


Thus,  i  n  order  to  i  nvesti  gatethesteady  state stabi  I  ity  theei  gen-val  ues  of  the  5x5  matrix  needs 
to  be  computed;  the  system  is  stable  iff  the  ei  gen-val  ues  have  negative  real  parts. 
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